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IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



In re application of DORNSTETTER, et al . 

U.S. patent application 
(U.S. National Phase of PCT/FR99/02089 

filed on September 2, 1999) 
(corr. to French applications 

No. 98 11095 and No. 98 12552) 



For : 11 DIGITAL EQUALISING METHOD AND RADIO COMMUNICATION 
RECEIVER IMPLEMENTING SAID METHOD " 



PRELIMINARY AMENDMENT 

To Honorable Commissioner of Patents and Trademarks 
Washington, D.C. 20231 



SIR: 

Prior to examination, please amend the above - 
identified application as follows: 



IN THE SPECIFICATION: 



Page 1, between lines 2 and 3, insert the heading 
BACKGROUND OF THE INVENTION 



Page 5, between lines 10 and 11, insert the heading 
— SUMMARY OF THE INVENTION 

Page 7, delete lines 29-32. 

Page 7, before line 33, insert the heading -- BRIEF 
DESCRIPTION OF THE DRAWINGS 

Page 8, between lines 3 and 4, insert the heading 
DESCRIPTION OF PREFERRED EMBODIMENTS 



IN THE CLAIMS: 

Please amend Claims 1-8 as follows: 

-- 1 (Amended). A digital equalization method for 

estimating discrete information symbols [ (d k ) of a 

transmitted signal] from digital samples [ (y k ) ] of a 

signal received over a transmission channel represented by 
a finite impulse response of W+l coefficients 
[ (r Q/ r 1# . . . , r w ) ] , W being an integer greater than 1, 

comprising the steps of: 

- determining [the] W roots [ (a x/ a 2 r . • . /O^) ] in the 

complex plane of the Z-transform [(R(Z))] of the 
impulse response; 

- distributing the W roots into a first set of W-p 

roots [ (a x , . . . , oc^p) ] and a second set of p roots 
[ (ct w . p+1 / . . . tQLff) ] / p being an integer greater than 0 
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and smaller than W, the roots of the second set 
being closer to the unit circle than those of the 
first set according to a determined distance 
criterion in the complex plane; 

- obtaining an intermediate signal C(Y')] by applying 
a first equalization method to the received signal 
[ (Y) ] based on a finite impulse response [whose] 
having a Z-transform [(R S (Z)), formed by] consisting 
of a polynomial of degree W-p in Z" 1 , [has] having 
roots [which are] egual to the W-p roots of the 
first set; and 

- obtaining estimations [ ( d k ) ] of the discrete 
information symbols [of the transmitted signal] by 
applying a second equalization method to the 
intermediate signal based on a finite impulse 
response [whose] having a Z-transform [(R X (Z)), 
formed by] consisting of a polynomial of degree p in 
Z' 1 , [has] having roots [which are] egual to the p 
roots of the second set. 

2 (Amended) . A method according to claim 1, wherein the 
first equalization method yields the intermediate signal 
in the form of a vector Y' of n+p samples [ (y' x # . . . #y' n+p ) ] 
obtained according to the relation : 

Y' = (A' H A' ) " x A' H Y 
where n is an integer representing a frame size, Y is a 
vector composed of n+W samples [ (y k ) ] of the received 
signal, and A' is a matrix with n+W rows and n+p columns 
having a Toeplitz structure formed from the coefficients 
[(s q )] of said polynomial of degree W-p in Z" 1 [(R S (Z))]. 
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3 (Amended) . A method according to claim 1 [or 2] , 
wherein the second equalization method comprises 
implementing a Viterbi algorithm, 

4 (Amended) . A method according to [any one of claims 1 
to 3] claim 1 , wherein the unit circle distance criterion, 
used to distribute the W roots o^,...,^ of the 
Z-transform [(R(Z))] of the channel impulse response into 
the first and second sets, is expressed as a distance 5 q of 



the form 8 q = 1 - 



q 



if 



< 1, and of the form 



5 a =1-1/ 



a,— 



if 



a q 



> 1, for l<q<W. 



5 (Amended) . A radio communications receiver comprising: 

- conversion means [(1,3,4)] to produce digital 
samples [ (y k ) ] from a radio signal received over a 
transmission channel represented by a finite impulse 
response of W+l coefficients [ (r 0 , r x , . . . , r w ) ] , W 
being an integer greater than 1; 

- means [(6)] for measuring the channel impulse 
response; 

- means for calculating [the] W roots [ (a x , a 2 , . . . , a^) ] 

in the complex plane of the Z-transform [(R(Z))] of 
the impulse response; 

- means for distributing the W roots into a first set 
of W-p roots [ (o^, . • • ,0^ ) ] and a second set of p 
roots [ (ct w . p+1 , . . . ,0^) ] , p being an integer greater 

than 0 and smaller than W, the roots of the second 
set being closer to the unit circle than those of 
the first set according to a determined distance 
criterion in the complex plane; 
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- a first equalization stage for producing an 
intermediate signal by applying a first equalization 
method to the received signal [ (y k ) ] based on a 

finite impulse response [whose] having- a Z- transform 
[(R S (Z)), forced by] consisting of a polynomial of 
degree W-p in Z" 1 , [has] having roots [which are] 
equal to the W-p roots of the first set; and 

- a second equalization stage for producing 
estimations [ ( d k ) ] of the discrete symbols of a 
signal carried on the channel by applying a second 
equalization method to the intermediate signal based 
on a finite impulse response [whose] having- a 
Z- trans form [(R X (Z)), formed by] consisting of a 
polynomial of degree p in Z" 1 , [has] having roots 
[which are] equal to the p roots of the second set, 

6 (Amended) . A receiver according to claim 5, wherein 
the first equalization stage is arranged to yield the 
intermediate signal in the form of a vector Y' of n+p 
samples [ (y' ± , . . wy' n+p ) ] obtained according to the 
relation: 

Y' = (A /H A' ) -1 A' H Y 
where n is an integer representing a frame size, Y is a 
vector composed of n+W samples [ (y k ) ] of the received 
signal, and A' is a matrix with n+W rows and n+p columns 
having a Toeplitz structure formed from the coefficients 
[{s q )] of said polynomial of degree W-p in Z" 1 [(R S (Z))]. 

7 (Amended) . A receiver according to claim 5 [or 6] , 
wherein the second equalization stage is arranged to 
implement a Viterbi algorithm. 
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8 (Amended) • A receiver according to [any one of claims 
5 to 7] claim 5 , wherein the means for distributing the W 
roots into the first and second sets make use of a unit 
circle distance criterion expressed as a distance 5 q of the 



form 
5 CT =1-1/ 



« q = 1- 



a, 



if 



a, 



< 1, 



and 



of 



the form 



a, 



if 



a 



q 



> 1, for l<q<w. 
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REMARKS 

Entry of the above amendment is respectfully requested. 

The amended claims 1-8 are substantially similar to original claims 1-8 as set out in the 
PCT International Application, with the exception that the new claims have been amended to 
omit multiple dependent claims and to otherwise conform to U.S. practice. 

Respectfully submitted, 

PIPER MARBURY RUDNICK & WOLFE 



Date: 





Michael L. Kenaga 
Reg. No. 34,639 



PIPER MARBURY RUDNICK & WOLFE 
P.O. Box 64807 
Chicago, IL 60664-0807 
(312) 368-4000 
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DIGITAL EQUALISING M ETHOD 
AND RADIO CO MMUNICATION RECEIVER IMPLEMENTING SAID METHOD 

The invention concerns the digital equalization of 
signals. It has important applications in the field of 
radio communications . 

The method is applied to the processing of a signal 
received from a transmitter over a channel between 
transmitter and receiver, whose response is known or has 
been estimated beforehand. One of the main issues is then 
that of the compromise between the performance of the 
equalizer and its complexity. 

A complete maximum likelihood estimation of the 
transmitted signal is possible, for example by employing 
the Viterbi algorithm (see 6.D. Forney Jr. : The Viterbi 
Algorithm, Proc. of the IEEE, Vol. 61, No. 3, March 1973, 
pages 268-278) . Nevertheless, as soon as the impulse 
response of the channels becomes long or the number of 
possible discrete values becomes large, the exponentially 
increasing complexity of these methods renders them 
impracticable . 

We consider the case of a radio communications 
channel used for the transmission of a signal composed of 
successive sequences or frames of n symbols, d k (l<k<n) . 
The symbols d k have discrete values : binary (±1) in the 
case of a BPSK (binary phase shift keying) modulation; 
quaternary (±l±j) in the case of a QPSK (quaternary phase 
shift keying) modulation. 

After baseband conversion, digitization and matched 
filtering, a vector Y of the received signal, 
corresponding to the symbols transmitted over the duration 
of a frame, is defined by the expression : 
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where W+l is the length in numbers of bits of the 
estimated impulse response of the channel, 

r = (r Q/ r x , . .., r w ) is the estimated impulse response of 
the channel, the r q being complex numbers such that r q =0 
for q<0 or q>W, y k is the k-th complex sample received 
with l<k<L=n+W, and Y N is a vector of dimension L composed 
of samples of additive noise y N/k . The estimated impulse 
response r takes into account the propagation channel, the 
signal shaping by the transmitter and the receiver 
filtering. 

The matrix A of dimension Lxn has a Toeplitz- type 
structure, meaning that if a i#j denotes the term in the 
i-th row and in the j-th column of the matrix A, then 
a i+l,j+l a a i,j for and l<i<n-l. The terms of the 

matrix A are given by : a X/J sO for l<j<n (A therefore has 

only zeros above its principal diagonal); a. -=0 for 

x , ± 

W+l<i<L (band-matrix structure); and a. 1= r. . for l<i<W+l. 

X , X X ~ _L 

The matrix equation (1) expresses the fact that 
the signal received Y is an observation, with an additive 
noise, of the convolution product between the channel 
impulse response and the transmitted symbols. This 
convolution product can also be expressed by its 
Z- trans form : 
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Y(Z) = R(Z) .D(Z) + Y N (Z) (2) 
where D(Z), Y(Z), R(Z) and Y N (Z) are the Z- trans forms of 
the transmitted symbols, of the received signal, of the 
impulse response and of the noise, respectively : 



D(Z) = J d k . Z~ k (3) 

L 

Y (2) = Zyk-Z" k (4) 
w 

R(Z) = £ r g . z " q (5) 

q=0 

A conventional solution to solve a system such as 
(1) is the so-called "zero forcing" method, by which we 
determine the vector £> ZF of n continuous components which 

• - • , _ , jl - j|2 

minimizes the quadratic error s = AD— Y . Subsequently, a 

discretization of the components of vector D ZF relating to 
each channel is performed, often through a channel 
decoder. The least-squares solution D ZF is given by : 

15 D ZF = (A H A) " X A H Y, where A H denotes the conjugate transpose of 
matrix A. We then are faced with the problem of inverting 
the Hermitian positive matrix A H A. The inversion can be 
effected by various classical algorithms, either directly 
(method of Gauss, Cholesky, etc.) or by an iterative 

20 technique (Gauss-Seidel algorithm, gradient algorithm, 
etc . ) . 

The estimation error D- D ZF is equal to (A H A) ~^*A H Y N , 
hence the solution includes noise with a variance: 

a 2 = e(|d-D zf | 2 ) = N 0 x Trace (a h a)~ 1 ] (6) 
25 where N Q is the noise power spectral density. It can be 



seen that noise enhancement occurs when the matrix A H A" is 
badly conditioned, i.e. when it has one or more 
eigenvalues close to zero. 

This noise enhancement is the main drawback of the 
conventional solution methods. In practice, the cases 
where the matrix A H A is badly conditioned are frequent, 
especially with multiple propagation paths. 

A relatively simple means to partly remedy this 
drawback is known, by accepting residual interference in 
the solution, i.e. by not adopting the optimum least- 
squares solution, but the solution : ^mmse" (^ H A+ N 0 ) -1 A H Y, 
where N Q denotes an estimation of the noise spectral 
density, that the receiver must then calculate. This is 
known as the MMSE (minimum mean square error) method and 
it allows the estimation variance to be reduced relative 
to the zero- forcing method, but introduces a bias. 

The zero forcing methods and the like amount to 
performing an inverse filtering on the signal received by 
a filter which models the transfer function 1/R(Z) 
calculated using a certain approximation (quadratic in the 
case of zero forcing) . Where one or more roots of the 
polynomial R(Z) (equation (5)) are situated on the unit 
circle, the theoretical inverse filter presents 
singularities rendering it impossible to estimate by a 
satisfactory approximation. In the case of the quadratic 
approximation, this corresponds to a divergence in the 
error variance a 2 when the matrix A H A has an eigenvalue 
equal to zero (relation (6)). 

This problem does not arise in methods such as the 
Viterbi algorithm which intrinsically take the discrete 
nature of the symbols into account, but which require a 
much higher calculating power for large- sized systems. 
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U.S. Patent 4,701,936 discloses a channel equalizer 
using an all -pass filter determined with reference to the 
Z- transform of the estimated channel impulse response. 

An object of the present invention is to propose an 
equalization method achieving a good compromise between 
the reliability of the estimations and the complexity of 
the equalizer. 

Another object is to produce an equalizer requiring 
a reasonable calculating power but capable of processing, 
with performance comparable to that of Viterbi equalizers, 
signals whose symbols have a relatively high number of 
states and/or signals carried on a channel with a 
relatively long impulse response. 

Accordingly, the invention proposes a digital 
equalization method for estimating discrete information 
symbols from digital samples of a signal received over a 
transmission channel represented by a finite impulse 
response of W+l coefficients, W being an integer greater 
than 1. This method comprises the steps of: 

- determining the W roots in the complex plane of the 
Z- transform of the impulse response; 

- distributing the W roots into a first set of W-p 
roots and a second set of p roots, p being an 
integer greater than 0 and smaller than W, the roots 
of the second set being closer to the unit circle 
than those of the first set according to a 
determined distance criterion in the complex plane; 

- obtaining an intermediate signal by applying a first 
equalization method to the received signal based on 
a finite impulse response whose Z-transform, formed 

by a polynomial of degree W-p in Z" 1 , has roots which 
are the W-p roots of the first set; and 

- obtaining estimations of the discrete information 
symbols by applying a second equalization method to 
the intermediate signal based on a finite 
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impulse response whose Z- transform, formed by a 
polynomial of degree p in Z" 1 , has roots which are 
the p roots of the second set. 

The "first equalisation method" will usually be 
5 chosen so as to treat the unknown symbols as continuous 
variables. This leads to a process similar to an inverse 
filtering whose transfer function would have a form 
approaching the expression 1/R S (Z), where R S (Z) denotes the 
polynomial in Z' 1 of order W-p having roots that are the 

10 W-p roots furthest from the unit circle in the complex 
plane. In particular, it could be of the "zero- forcing" 
type. This process generates only a small enhancement of 
the noise, since the roots of the associated transfer 
function in Z are relatively far from the unit circle. 

15 For the p roots which are closest to the unit 

circle, we adopt measures to avoid or to limit the 
incidence of the noise enhancement problem. A MMSE or 
similar method can be chosen as the "second equalization 
method". However, this second method will advantageous 

2 0 take into account the discrete nature of the unknown 

symbols. In particular, it may rely on a trellis 
algorithm, such as the Viterbi algorithm, whose 
implementation is common in channel equalizers where the 
system size is not too big. 
25 The implementation of the second equalization method 

is generally more complex than the first. In each 
particular case, the choice of the number p allows the 
best compromise between the reliability of the 
estimations, which favors higher values of p, and the 

3 0 complexity of the equalizer, which favors lower values of 

P- 

Another aspect of the present invention relates to a 
radio communications receiver comprising : 



- conversion means to produce digital samples from a 
radio signal received over a transmission channel 
represented by a finite impulse response of w+1 
coefficients, W being an integer greater than 1; 

- means for measuring the channel impulse response; 

- means for calculating the W roots in the complex 
plane of the Z- transform of the impulse response; 

- means for distributing the W roots into a first set 
of W-p roots and a second set of p roots, p being an 
integer greater than 0 and smaller than W, the roots 
of the second set being closer to the unit circle 
than those of the first set according to a 
determined distance criterion in the complex plane; 

- a first equalization stage for producing an 
intermediate signal by applying a first equalization 
method to the received signal based on a finite 
impulse response whose Z transform, formed by a 
polynomial of degree W-p in Z' 1 , has roots which are 
the W-p roots of the first set; and 

- a second equalization stage for producing 
estimations of the discrete symbols of a signal 
carried on the channel by applying a second 
equalization method to the intermediate signal based 
on a finite impulse response whose Z transform 
(R X (Z)), formed by a polynomial of degree p in Z" 1 , 
has roots which are the p roots of the second set. 
Other features and advantages of the present 

invention will become apparent in the following 
description, with reference to the attached figures, of 
non- limiting embodiments, in which : 

- Figure 1 is a block diagram of an example of a radio 
communications receiver according to the invention; 
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- Figure 2 is a flow chart showing an embodiment of 
the method according to the invention; and 

- Figure 3 is a diagram illustrating the performances 
of the method. 

5 The receiver represented in Fig. 1 comprises a radio 

stage 1 which receives the signal picked up by antenna 2 
and converts it to the baseband. The baseband signal is 
digitized by the analog- to-digital converter (ADC) 3 and 
fed to the reception filter 4. Filter 4 provides matched 

10 filtering with respect to the signal shaping by the 
transmitter. It delivers a digital signal at a rate of one 
complex sample per transmitted symbol. 

The digital signal is fed to a demodulator 
comprising, on the one hand, a synchronization and channel 

15 estimation module 6 and, on the other hand, an equalizer 
7. 

The synchronization and channel estimation are, for 
example, performed conventionally by means of a 
synchronization sequence included by the transmitter in 
20 each signal frame. Detection of this sequence, known to 
the receiver, allows, on the one hand, the receiver to be 
synchronized relative to the temporal structure of the 
transmitted frames, and on the other to estimate the 
impulse response r = (r Q , r 1# . . . , r w ) of the channel carrying 

25 the transmitted frames. The impulse response calculated by 
module 6 is fed to the equalizer 7. 

Equalizer 7 functions for example according to the 
flow chart shown in Fig. 2 to process each frame of the 
received signal, which are presented in the form of a 



3 0 vector Y= 



, with L=n+W, adopting the previous notation. 



The channel estimation module 6 having supplied the 



W+l complex coefficients r q of the estimated channel 

impulse response, the first step 10 consists of a search 
of the W roots of the Z- transform of the impulse response, 
given by equation (5) . Various conventional methods for 
5 the determination of complex roots of a polynomial can be 
employed at step 10* To this end, reference can be made to 
the work of E. DURAND : ^Numerical Solutions of Algebraic 
Equations, Volume I, Equations of the Type F(x)=0", 
Editions Masson, 1960. 
10 The W roots thus found a 1 ,a 2 , . . . ,0^ are then ordered 

to allow distribution into two sets, one containing the 
W-p roots furthest from the unit circle, and the other the 
p roots closest to the unit circle. 

To this end, a distance 5 is calculated at step 11 

15 for each of the roots a q (l<q<W) . This distance can be 



conveniently obtained as follows 



5 q = i 



1 

1-1/ 



a q 



S 2. 

si 



a q 



<1 

(7) 

>1 



In step 12, the roots a q of the transfer function 
R(Z) are sorted in order of decreasing distance: 
2 0 8 X > 8 2 > . • • > 8 W . Then, the first W-p roots a 1# . . . /Ow_p' 
which are the furthest from the unit circle, are separated 
from the remaining roots ot w _ p+1 , . . . ,0^. 

In step 13, the equalizer 7 develops a polynomial in 
Z" 1 defined by : 

w-p w-p 

25 R S (Z) = J! i 1 - a g- 2_1 ) = X s q- (8) 

q=l q = 0 

This allows the coefficients s g of the transfer 

function R S (Z) associated with the impulse response 
£5 = (s Q , s x , . . . , s w _ p ) of a virtual channel to be determined, 
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10 



15 



20 



which corresponds to the estimated transmission channel 
with the contributions closest to the singularity zones 
removed. 

We then proceed to a first equalization 14 which is 
effectively an inverse filtering approaching the transfer 
function 1/R S (Z). Several implementations can be used to 
perform this inverse filtering. In particular, a zero- 
forcing equalization can be performed, as previously 
discussed. Regarding these methods, reference can be made 
to the work of J.G. Proakis : "Digital Communications'', 
McGraw-Hill, 2 nd Edition, 1989. 

The inverse filtering 14 produces an intermediate 
signal in the form of a vector Y' of L'=n+p samples 
Y f X f # * • 'y'li* In fc ke case of the zero-forcing method, the 
vector Y r is obtained from the matrix equation : 

Y' = (A ,H A' ) - 1 A' H Y (9) 
In expression (9) , A' denotes a matrix having a 
Toeplitz structure with n+W rows and n+p columns, 

generated from the coefficients s q of the polynomial R S (Z): 



*1 



0 
0 



0 
0 



The 



>w-p 

0 0 s w _j 
0 0 

sorting of the roots 



0s w- P y 



(10) 



results in the 



25 



eigenvalues of matrix A' H A' being relatively far removed 
from 0 . 

Alternatively, the inverse filtering can be achieved 
by cascading W-p filtering cells, each corresponding to 



11 



the inverse of a transfer function R| (Z) = 1 - c^Z" 1 , for 
l<q^W-p. if | a J = i # the inverse filter of R* (Z> is 

unfeasible. If | a q | < 1, we can develop 1/R| (z) in the 
form : 

~^ = 1+ a q . Z" 1 + a 2 . z~ 2 + . . . + a*. Z" m + . . . (ii) 

Development (11) i s causal and stable since the 

convergence domain contains the unit circle. The inverse 

filter cell can thus be generated in either a transverse 
or a recursive form. 

If |a q | > 1, then 1/r| (z) can be developed in the 

form : 

= -a, 1 , z . (l + c^ 1 . z + a~ 2 . z 2 + . . . + a~ m . z m + . . .) (i 2 ) 

Development (12) is anti-causal and stable. To perform the 
inverse filtering cell, the development (12) is truncated 
and an transverse filter implementation is adopted. The 
anti -causality results in a delay corresponding to the 
length of the retained response. 

We note that developments (11) and (12) justify the 
unit circle distance criterion 8 q used according to 
relation (7) . 

In step 15, the equalizer 7 develops a polynomial of 
degree p in Z" 1 , whose roots correspond to the p roots of 
R(Z) closest to the unit circle, such that 
R(Z) = R s (z) .R I (Z) : 



R X (Z) = r 0 . II (l " cc g . Z" 1 ) = X V Z" q (13) 
q=W-p+l q=0 



The complex coefficients t q define the impulse 
response of another virtual transmission channel, whose 
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equalization by a zero- forcing method or the like would 
pose problems of noise enhancement. 

The intermediate signal Y' is then subjected to an 
equalization according to another method, based on the 
impulse response t = <t 0 , t x , . . . , t p ) . This second 
equalization 16 can be conveniently performed by means of 
a Viterbi trellis (see above-cited article from G.D. 
Forney Jr., or the above-cited work of J.G. Proakis) . 

The second equalization stage 16 generates 
estimations d k of the symbols of the frame (l<k<n) . These 
estimations d k constituting the output of equalizer 7 can 
be input to a deinterleaving module 8 then to a channel 
decoder 9 which detects and/or corrects possible 
transmission errors. 

Figure 3 illustrates the performance of the method 
in the case of the transmission of a signal frame 
according to the GSM European mobile telephone system, 
replacing the binary modulation of the GMSK type by an 
eight state phase modulation (8-PSK modulation) . The 
impulse response was truncated at five bit time periods 
(W=4) . Figure 3 shows the dependence of the bit error rate 
(BBR) , expressed in %, on the ratio E b /N Q between the 
energy-per-bit and the spectral density of the noise, 
expressed in decibels (dB) . The BER is that observed in 
the symbol estimations after deinterleaving and channel 
decoding carried out by a method conforming to those 
employed in the GSM standard. Curve I shows the results 
obtained by the pure zero- forcing method, i.e. in the 
limit case where p=0. Curve II shows the theoretical 
result that would be obtained by equalizing the channel 
purely with the Viterbi algorithm (limit case where p=w) . 
In practice, the corresponding trellis should have 8 4 =4096 
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states, meaning that the Viterbi equalizer could not be 
realized using current technology. The difference between 
curves I and II illustrates the superiority of the Viterbi 
algorithm which produces maximum likelihood estimations. 
5 Curves III and IV show the results obtained by the 

method according to the invention, in the cases where p=l 
and p=2, respectively. The quite significant improvement 
already achieved for the value p=l with respect to pure 
zero forcing can be seen. 

10 By way of an example , equalization of a GSM signal 

frame by a pure Viterbi algorithm, under the conditions of 
Figure 3, would require of the order of 8.45 million 
floating-point operations, i.e. around 1.83 Gflops, while 
the implementation of the present invention under the same 

15 conditions requires of the order of 19,000 floating-point 
operations («4.2 Mflops) in the case of p=l, including 

searching the roots of R(Z) and 1/R S (Z) inverse filtering 
by the zero- forcing method. The number increases to the 
order of 129,000 operations («28 Mflops) in the case of 
20 p=2, which remains compatible with the capacity of 
currently available digital signal processors (DSP) . 
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CLAIMS 

1. A digital equalization method for estimating 

discrete information symbols (d k ) from digital samples (y k ) 

of a signal received over a transmission channel 
represented by a finite impulse response of W+l 
coefficients (r Q , r^, • . . , r w ) , W being an integer greater 

than 1, comprising the step of determining the W roots 
(0^, a 2 , . . . , c^) in the complex plane of the Z-transform 

(R(Z)) of the impulse response, characterized it further 
comprises the steps of: 

- distributing the W roots into a first set of W-p 
roots (a lf • . ♦ , oc^p) and a second set of p roots 

(a w _ p+1/ . ♦ . ,0^) / p being an integer greater than 0 

and smaller than W, the roots of the second set 
being closer to the unit circle than those of the 
first set according to a determined distance 
criterion in the complex plane; 

- obtaining an intermediate signal (Y') by applying a 
first equalization method to the received signal (Y) 
based on a finite impulse response whose Z-transform 

(R S (Z)), formed by a polynomial of degree W-p in Z" 1 , 
has roots which are the W-p roots of the first set; 
and 

- obtaining estimations ( d^ ) of the discrete 
information symbols by applying a second 
equalization method to the intermediate signal based 
on a finite impulse response whose Z-transform 

(R I (Z)) / formed by a polynomial of degree p in Z" 1 , 
has roots which are the p roots of the second set. 
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2, A method according to claim 1, wherein the first 
equalization method yields the intermediate signal in the 
form of a vector Y' of n+p samples (y ' x , . . * ,y ' n+p ) obtained 

according to the relation : 

Y' = (A' H A')" 1 A' H Y 
where n is an integer representing a frame size, Y is a 
vector composed of n+W samples (y k ) of the received 

signal, and A' is a matrix with n+W rows and n+p columns 
having a Toeplitz structure formed from the coefficients 
(s q ) of said polynomial of degree W-p in Z" 1 (R S (Z)). 

3, A method according to claim 1 or 2, wherein the 
second equalization method comprises implementing a 
Viterbi algorithm, 

4, A method according to any one of claims 1 to 3, 
wherein the unit circle distance criterion, used to 
distribute the W roots o^,...,^ of the Z-transform (R(Z)) 
of the channel impulse response into the first and second 
sets, is expressed as a distance 5 q of the form 



8 q - 1 - 



1 



if 



a q 



if 

> 1, for l<q<W 



^ 1, and of the form 6 q =1-1/ 



A radio communications receiver comprising : 

- conversion means (1,3,4) to produce digital samples 
(y k ) from a radio signal received over a 

transmission channel represented by a finite impulse 
response of W+l coefficients (rg,^, , . . ,r w ) , W being 

an integer greater than 1; 

- means (6) for measuring the channel impulse 
response; 
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- means for calculating the W roots (o^, cc 2 , . . . , o^) in 

the complex plane of the Z-transform (R(Z)) of the 
impulse response ; 

- means for distributing the W roots into a first set 
of W-p roots (ot 1/ ...,a w _ ) and a second set of p 
roots (a w-p+1 , . . . , c^) , p being an integer greater 
than 0 and smaller than W, the roots of the second 
set being closer to the unit circle than those of 
the first set according to a determined distance 
criterion in the complex plane; 

- a first equalization stage for producing an 
intermediate signal by applying a first equalization 
method to the received signal (y k ) based on a finite 

impulse response whose Z transform (R S (Z)), formed 
by a polynomial of degree W-p in Z" 1 , has roots which 
are the W-p roots of the first set; and 

- a second equalization stage for producing 
estimations ( d k ) of the discrete symbols of a signal 
carried on the channel by applying a second 
equalization method to the intermediate signal based 
on a finite impulse response whose Z transform 
(R X (Z)), formed by a polynomial of degree p in Z" 1 , 
has roots which are the p roots of the second set. 

6. A receiver according to claim 5, wherein the first 

equalization stage is arranged to yield the intermediate 
signal in the form of a vector Y' of n+p samples 
<Y'l* • • • 'Y'n+p) obtained according to the relation : 

Y' = {A' H A') - 1 A /H Y 
where n is an integer representing a frame size, Y is a 
vector composed of n+W samples (y k ) of the received 
signal, and A' is a matrix with n+W rows and n+p columns 
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having- a Toeplitz structure formed from the coefficients 
£s q ) of said polynomial of degree W-p in Z' 1 (R s (z)), 

7. A receiver according to claim 5 or 6, wherein the 
second equalization stage is arranged to implement a 
Viterbi algorithm. 

8, A receiver according to any one of claiuas 5 to 7, 
wherein the means for distributing the W roots into the 
first and second sets make use of a unit circle distance 
criterion expressed as a distance 5 q of the form 

5q =* 1 - |a q | if ja^l ^ 1* and of the form S q = 1 - l/Jctq) 

if |a q | > 1, for l<q£W. 
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